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Abstract 

It is proved that irrespective of compressibility tokamak steady states with 
purely poloidal mass flow can not exist in the framework of either magneto- 
hydrodynamics (MHD) or Hall MHD models. Non-existence persists within 
single fluid plasma models with pressure anisotropy and incompressible flows. 
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I. Introduction 



Motivation of the present study was a proof that ideal magnetohydrody- 
namics (MHD) steady states of magnetically confined plasmas with purely 
poloidal incompressible flows and magnetic fields having toroidal and poloidal 
components can not exist [H12J- The inconsistency relates to the toroidicity 
because in cylindrical geometry equilibria with purely poloidal (azimuthal) 
flows exist without restriction on the direction of the magnetic field [3] . The 
aforementioned statement of non-existence includes tokamak equilibria; only 
field reversed configurations with purely poloidal incompressible flows which 
by definition have only poloidal magnetic fields are possible. Also, magnetic 
dipolar steady states with purely poloidal compressible flows parallel to the 
magnetic field were studied in Ref. [I]. In tokamaks, however, poloidal 
sheared flows play an important role in the transition from low to high con- 
finement mode (L-H transition). In fact, in many experiments no toroidal 
rotation in connection with the transition is reported (see for example Ref. 
[5]). It is therefore interesting to theoretically examine wether alternative or 
additional physical input to incompressibility can remove the incompatibility. 

The present study aims at examining whether (i) compressibility (ii) two 
fluid effects and (iii) pressure anisotropy can give rise to the existence of 
tokamak steady states with purely poloidal flows. It turns out that in cases 
(i) and (ii) as well as in case (iii) for incompressible flows the non-existence 
statement keeps hold. It should be noted that there is a number of papers 
on equilibria with flow within the framework of MHD [2]- [I], [S]-[IZ], the 
two fluid model [?J HU [17], [18]-[2T] and for anisotropic pressure [22]-[25]. 
Certain of the derivations and equations therein are related to the present 
work. We will prefer, however, to present the study in a self contained way 
because explicit reference to the content of the aforementioned papers for 
the three models employed would make a common notation difficult and the 
manuscript inconveniently readable. Compressible axisymmetric MHD equi- 
libria is the subject of Sec. II. In Sec. Ill two fluid effects are examined in the 
framework of Hall MHD model with electron temperatures uniform on mag- 
netic surfaces and either ion temperatures uniform thereon or incompressible 
ion flows. Pressure anisotropy for incompressible flows is considered in Sec. 
IV. The conclusions are summarized in Sec. V. 
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II. Compressibility 



The ideal MHD equilibrium states of a magnetically confined plasma are 
governed in standard notation and convenient units by the following set of 
equations: 

V ■ (pv) = 0, (1) 
p (vV)v = jxB-VP, (2) 

V x E = 0, (3) 
VxB=j, (4) 

V • B = 0, (5) 

E + vxB = 0. (6) 

An energy equation or equation of state needed to close the set of Eqs. ([!])- 
(EJ) is not adopted from the beginning; it will be specified when necessary 
later. On account of axisymmetry and Ampere's law (00) the divergence-free 
fields, i.e. the magnetic field B, the current density j, and the momentum of 
the ion fluid element pv can be expressed in terms of scalar functions ?p(R, z), 
I{R,z), F(R,z) and Q(R, z) as 

B = I(R, z)V0 + x 'Vip(R, z), (7) 

j = A*^V0 — V0 x V/(-R, z), (8) 

pv = Q{R,z)V<j) + V0x VF(R,z). (9) 

Here, (R, z, 0) are cylindrical coordinates with z corresponding to the axis of 
symmetry; the functions ip and F label the magnetic and velocity surfaces, 
respectively; A* = i? 2 V ■ (V/i? 2 ). Expressing the electric field in terms of 
the electrostatic potential, E = — V$, the component of Ohms law along 
V0 implies that function F is uniform on magnetic surfaces: 

F = F(i/j). (10) 

The electrostatic potential is also a surface quantity 

$ = $(V0, (ii) 
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as it follows by projecting along B. Two other integrals are identified 
in terms of surface quantities by the components of (jSD along V^>, and fl2]) 
along B, respectively: 

{IF' — 6) = (12) 




/ 1-^- )+R 2 F'^' = X(ip). (13) 



Eqs. ( TT2l and ({TBI can be solved for the functions I and associated with 
the toroidal components of B and v to yield 

X - R 2 & 

= (14) 

r,, , ™ XF'-R 2 p& . 
e <**> = I -Ml ■ < 15) 

Here M 2 is the square of the Mach function of the poloidal velocity with 
respect to the Alfven velocity: 



M 
p 



'\2 



VpAj P 



In the present study we will consider both compressible and incompressible 
flows. For typical tokamak poloidal velocities and temperatures (maxt> p pa 
10 4 m/sec and max kTi pa 10 keV) the incompressibility condition 



Vtht 



« 1, (16) 



where Vth t is the ion thermal velocity, is satisfied and therefore incompressible 
flows are of relevance. Also, it may be noted that incompressibility is a good 
approximation for small flows lying within the first elliptic regime of the 
equilibrium differential equations. In MHD this regime is 

< M p < /3 (17) 

(see for example Ref. [9]), where the maximum value of j3, defined as the ratio 
of the plasma pressure to the total (plasma and magnetic field) pressure, is 
about 0.35 (taking into account the particularly high values of j3 which have 
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been obtained in spherical tokamaks). Note that (1171) sets nearly the same 
bounds on v p as (TTSj) as can be seen by using the definitions of M p and j3. 

On the basis of fTH| and ffl5l) let us first recover the non-existence of equi- 
libria with purely poloidal incompressible flows and magnetic fields having 
poloidal and toroidal components. On account of incompressibility, V-v = 0, 
([!]) implies that p = p(ip). Then, for 9 = 0, f[T51) can not be satisfied, be- 
cause in addition to surface quantities it contains R explicitly, unless either 
F' = $' = or X = $' = 0. The former case corresponds to a static 
equilibrium (v = 0). The latter implies that v is parallel to B and there- 
fore = as it also follows from (fl4|) . Therefore, tokamak equilibria with 
purely poloidal incompressible flows are not possible. Note that if, inversely, 
the magnetic field is purely poloidal (J = 0) inspection of Eqs. (TT^|) and (FTSl) 
implies that the flow is either purely poloidal or purely toroidal. Therefore 
in field reversed configurations coexistence of toroidal and poloidal incom- 
pressible flows is not admitted. 

In the case of compressibility the density can vary on magnetic surfaces. 
For purely poloidal flows, however, Eqs. (I14p and (fT51) imply / = X(i/j) and 
restrict p to be of the form 

XF' , . 

" = vw (18) 

Evaluation of Vp on magnetic axis, on which V^> = 0, by (118j) yields 

V Po = -2Rl^e R . (19) 

Therefore, 

Vp + 0, (20) 

unless X Q = or Fq = which would imply p = by (|T8|) . To proceed 
further we need an energy equation or equation of state. Since thermal 
conductivity along B is very large in high temperature plasmas, isothermal 
magnetic surfaces is an appropriate equation of state for tokamaks. If T = 
T(ip), using the ideal gas law P = apT the component of ([2]) along B yields 
for 6 = 

aTlnp = 0. (21) 




or 

y + aTlnp = H(i;). (22) 
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Evaluation of Vp on axis from (1221) after multiplying it by p and using Q 
for v, G = and (I18p for p yields 



[# -aT (l + lnp )]Vpo = 0. 



(23) 



If Hq = aTo(l+ln po) it follows from fl22|) that Tq = 0. Therefore, (1231) implies 
Vpo = which contradicts (1201) . As can be shown along the same lines 
the non unique definition of Vp on axis persists if alternative equations of 
state are adopted such as isentropic magnetic surfaces or barotropic plasmas 
[P = P(p)]. Thus, we can conclude that the non existence of tokamak 
equilibria with purely poloidal flows is extended to the compressible regime. 

The above local proof of non existence can be extended near the magnetic 
axis. Indeed, on account of (ITS]) . Eq. (1221) takes the form 



(RFJ 
2G 2 



l\2 



'df\ (Of 
dR + [ ~dz~ . 



(24) 



where G(ip) = XF'/&. The /^-derivative of ( 1241) on axis yields 



(jtq 



dip d 2 tp 



dR OR 2 



G 



(25) 



Note that the term involving d 2 ip/dRdz\ is not included in (125]) because this 
derivative vanishes as it follows from a Mercier expansion based consideration 
of the generalized Grad-Shafranov equation [Eq. ( |27j) below] around the 
magnetic axis. Since the RHS of (12"51) is finite, d 2 ip/dR 2 \ must tend to 
infinity. Consequently, the validity of the statement follows on the basis of a 
Mercier expansion of any solution of (pH]) near the magnetic axis: 



ip(x,y) =ip + 



1 d 2 i) 



2 8R 2 



2q, Ro 



+ 2 dz 2 



y 2 + 



(26) 



Here, (x, y) are Cartesian coordinates defined by R = Ro + x and z = z + y 
and ipo refers to the magnetic axis. Also, the current density on axis becomes 
singular. Furthermore, it may be noted that projection of ([2]) onto yields 
the generalized Grad-Shafranov equation 



- M l) 



2\ W 



R 2 



+ 



F'F" |V^| 2 IX' 



R 2 



R 2 



+ pH' + ap{l-\np)T = 0. (27) 
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Therefore, ip should satisfy the two differential equations (}2"Tj) and (J27J) con- 
taining the surface quantities T(ijj), H(ijj), F(ijj) and X(ip). As will be shown 
in the next section for incompressible flows the solution of these equations is 
irrelevant to tokamaks. 

III. Two fluid effects 

We will employ the Hall MHD, a simple two- fluid model in the approxi- 
mation of very small electron mass. Consequently, the electron momentum 
equation can be put in the Ohm's law form 

E + vxB = -(jxB - VP e ) , (28) 
P 

where v is the ion fluid element velocity, p = nM and h = M/e. The right 
hand side (RHS) of ( |28|) contains the Hall and electron pressure gradient 
terms. Eq. ( |28|) replaces ([6]) which is formally recovered in the limit of 
h — > 0. The other equations of the model are identical with (CQ), (T5}|5]) on 
the understanding that v = Vj. The momentum equation ([2]) is derived by 
a superposition of the electron and ion momentum equations neglecting the 
electron convective velocity term because of the very small electron mass. 
As in Sec. II certain integrals will be identified in the form of conserved 
quantities on magnetic surfaces by projecting (J6J) and (J2J) onto V0, B and 
V-^. The axisymmetric representations ©-Q of the divergence free fields 
will be employed in the derivations to follow which up to Eq. f|33|) below hold 
for velocities of arbitrary direction. Assuming that the electron temperature 
is uniform on magnetic surfaces and using E = — V$ and P e = apT e (ip) the 
component of §5§ along B yields 

B ■ V ($ - haT e log p) = (29) 
$ - haTJogp = E(ip). (30) 



or 



Eq. (1281 then becomes 

v x B = -j x B + - ha— - (1 - logp) V0. (31) 

p dip dip 

Projecting (fT3"|) along yields another integral: 

F + hI = fty). (32) 
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Eq. f)32p implies that in general the velocity surfaces depart from magnetic 
surfaces. The fact that v and B share the same surfaces in MHD [Eq. fllOp ] 
is recovered for h — > 0. Furthermore, the component of (13~T1) along leads 
to the elliptic differential equation 



hA*i/j = pR 2 



For h — > (1331 reduces to the algebraic MHD equation ([12]) . Purely poloidal 
ion flows will be further considered either incompressible or compressible on 
an individual basis as follows. 

Incompressible ion flow 

It is convenient to introduce the generalized vorticity 

ft = B + hV x v. (34) 

For purely poloidal flows the il-surfaces coincide with the magnetic surfaces. 
Since f2 is divergence free it can be expressed by 

Cl = N(R,z)V0 + V<px Vip. (35) 

Using (J31D, P = P e + Pi, E = -V$ and the identity 

v 2 

(v ■ V) v = V v x V x v, 

Eq. §2§ can be cast in the form 

pVH> = pv x Q - hVP u (36) 

where 

v 2 

W = h— + $. (37) 

The component of ( 1361) along V0 yields 

F = F(ip). (38) 

Eq. fl38l) implies that v lies on magnetic surfaces and therefore p = p(ip) 
because of incompressibility. Note that the electron fluid element velocity lies 
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on magnetic surfaces too whatever is the direction of v e as it follows from the 
electron momentum equation with the (v e • V) v e term being neglected. From 
([30]) and ((32} then it follows $ = $(^) and / = I (iff). Also, substituting 
the expressions for B and v from ([7]) and (Q into fl35|) leads to the following 
expression for N: 

fF'\ F' 
N = I + hV — ■ VV> + h—A*ip. (39) 

VP/ P 



Projecting (1361) onto fi and V?/>, respectively, and using ( 1391) for A" furnishes 

/iP< = hPuW ~ pW, (40) 

kMjWAV + ^ W = if (/^ - - (41) 

Here, the surface quantity -Pj S (V0 for vanishing flow coincides with the ion 
pressure and M£ = (F') 2 /p. Eqs. ( )4T!) and (|33|) for 6 = can be cast in the 
forms 

A*ii> = -fty)-R 2 gU>), (42) 

|V^| 2 = 2[2(^) + i? 2 j'W], (43) 

where f,g,i and j are known functions of p, F',T e , Pi S $ and /. The forms of 
Eqs. ( 142 p and (I43p indicate that the magnetic surfaces are identical in shape 
with those of the Palumbo solution [26] which, however, can not describe 
tokamak equilibria (iff contours of this solution are provided in Figures 1 and 
2 of Ref. [27]). Note that the pressure is not uniform on magnetic surfaces 
and therefore the equilibrium is not isodynamic. Only under the additional 
assumption P = P(i[f), (1401) implies B 2 = B 2 (ip). Recapitulating, two fluid 
effects in the frame of Hall-MHD model result in non pertinent to tokamaks 
isodynamic like equilibria with purely poloidal incompressible flows. 



Compressible ion flow 



In connection with tokamak equilibria we will assume that additionally 
to the electron temperature the ion temperature is uniform on magnetic sur- 
faces. Some of the derivations in the previous part of this section, i.e. Eqs. 
(J5D]) . f[3"2"j) . (13"3"j) and (13"6"1) . remain valid. Also, for purely poloidal flows irre- 
spective of compressibility the velocity surfaces coincide with the magnetic 
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surfaces [F = F(ip)] and therefore / = I(ip) by (13"2|) . Using Pi = apT^ip) the 
component of (13"6"j) along B yields 



H> + /iaT;lnp = A(^). (44) 

From (1301 . ( 1371) and (1441) eliminating the functions $ and W we obtain for 
= 

fr l ; ' 1 = AW - Sty) - /ia(T e + Tj) In p. (45) 

Acting the gradient operator on (1431) and evaluating the resulting equation 
on magnetic axis leads to 

Tm + T e0 

Vp = 

Po 

and therefore Vp = 0. On the other side, on account of (J3"2j) and / = /(^>), 
(l3"3"j) for 6 = becomes 

hA*ij> = pR 2 [E' - haT' e (1 - logp)] - (F' + hi') I. (46) 

Evaluating the gradient of (1461) on magnetic axis leads to 

R 2 (H' + haT' e0 lnp ) Vp = -2i? Po2[ ) e /? +/i [VA^ + 2aR T^ oPo (1 - lnp ) e R ] 

(47) 

A prerequisite that (14T|) is compatible with Vp = is that the RHS of 
this relation vanishes. The RHS of (|47j) . however, consists of the first large 
MHD-like term and the second term on the order of h which for tokamaks 
is small. Therefore, these terms can not cancel each other. Note that the 
first term involving the product po^g remains always finite because if Sq 
becomes very small, this in the MHD-like limit of h — > would imply very 
large densities by (|T8l) which keeps valid in this limit. Therefore, it follows 
from (j47l) that Vpo 7^ 0. It may be noted that the assumption of weak 
enough two-fluid equilibrium effects leading to non- unique definition of Vpo 
on axis is fulfilled in tokamaks particularly in the major interior part of the 
plasma. The inconsistency in connection with Vpo persists if alternative ion 
equations of state are adopted such as barotropic [Pi = Pi(p)} or isentropic 
ion velocity surfaces. Therefore, as in MHD, tokamak equilibria with purely 
poloidal compressible flows in the framework of Hall-MHD model can not 
exist. 
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IV. Pressure anisotropy 



We will employ an one fluid model with the momentum equation ([21 
replaced by 

p(v ■ V)v = j x B — V • "P, (48) 

where the tensor 

V = P\\1 + aBB (49) 
is associated with the pressure anisotropy, a measure of which is the quantity 

<'= JL ff-- ( 5 °) 

The other equations (CQ) and (G2)-([6]) remain unchanged. Double adiabatic 
equilibria can be obtained in the framework of the Chew Goldberger Low 
equations of state [30]: P\\ = S\\(ip)p 3 / B 2 and Pj_ = S±(ip)pB where S\\ and 
S± are arbitrary surface quantities. Also, adiabatic MHD equilibria with 
isotropic pressure can be recovered by setting P\\ = P± = S(ijj)p y , where S 
is the specific entropy. Here we will consider incompressible flows without 
other specification of P\\ and P±. The integrals (I10l) - fll2l) stemming from the 
Ohm's law remain valid. The component of (I48p along B yields 

B ■ V \l (l - Ml - a) + R 2 F'&] = (51) 



or 

/ h _ (Q! _ ^ + R 2 F'& = X(iP). (52) 

For isotropic pressure [(a = 0)], fl52|) reduces to ([TBI) . As in the isotropic case 
( fT2l) and ( |52|) can be solved for / and to yield 

X-R 2 ® . . 

(53) 



1 - M 2 - a 



For purely poloidal flows ( [531) and ( |54|) imply 



P = RWTl - a) (55) 
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and 

/ = jiL. (56) 

If a = a (ip) it follows from fl55l) that equilibria with purely poloidal 
incompressible flows are not possible (when ^ and B p ^ 0) because 
then these flows are not compatible with p = p(ip) . Static equilibria with 
a = a(i/)) were investigated in Refs. [2H1 EH]- Compatibility of (1531) with 
p = p("0) restricts a to be of the form 

where /(^>) is an arbitrary surface quantity. Eq. fl56|) then implies / = 
R 2 X (tp) / f (tp) . Apparently, possible equilibria of this kind are irrelevant 
to tokamaks in which the toroidal magnetic field, = I/R, must have a 
vacuum component proportional to Rr 1 . 



V. Conclusions 



We have shown that the non-existence of ideal MHD tokamak equilib- 
ria with purely poloidal incompressible flow can be extended to the cases 
of (i) compressible MHD flows (ii) Hall-MHD incompressible and compress- 
ible flows and (iii) one fluid equilibria with pressure anisotropy and incom- 
pressible flows. The non-existence relates to the toroidicity. Specifically, for 
incompressible flows an inconsistency appears in terms of relations which 
in addition to surface quantities have an explicit dependence on R or only 
iso dynamic-like equilibria are possible and for compressible flows the density 
gradient does not have unique definition on the magnetic axis. Also, for 
MHD compressible flows the proof can be extended near the magnetic axis. 

The existence of this kind of equilibria on account of additional physi- 
cal input, e.g. in plasmas with pressure anisotropy and compressible flows, 
or other models remains an open question. Experimental clarification on 
whether an even small toroidal velocity is necessary for the L-H transition to 
occur would be helpful. If not necessary, the existence of tokamak equilibria 
with purely poloidal flows from the theoretical point of view should remain 
an even more challenging problem. 
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